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Abstract. We find an upper bound for geodesic distances associated to 
monotone Riemannian metrics on positive definite matrices and density ma- 
J> ■ trices. 

o 

The notion and importance of Fisher information is well established 
in statistics and probability theory. As a measure of distinguishabil- 
ity of probability densities, the Fisher information was used by Rao 
to define a Riemannian metric on probability spaces. On the simplex 
of probability vectors V n = {p = (pi, . . .,p n ), Y^iPi = l >Pi > 0,i = 
1, . . . , n}, this is the unique metric contracting under markovian map- 
pings, by the Chentsov uniqueness theorem. On V n , the Fisher metric 
is 

y) = S ^ j Pi 1 x l y i x, y G T p V n 

i 

The geometry of V n with this metric is quite simple. By 

(1) P^ 2(VPT, ...,^p~n), 

it is isometric with an open subset in the sphere of radius 2 in M. n , 
The metric can be extended to the set A4 n — {p — (pi, . . . ,p n ), Pi > 0} 
of all finite (strictly positive) measures on the set {l,...,n}. Using 
the isometry f^Q) and elementary geometry in R n , we may compute the 
geodesic distance for the Fisher metric in V n and M. n : 

D(p, q) = 2 arccos(J^ ^/pl^/qT) p,q G V n 

i 

(the Bhattacharya distance) and 

d{p, q) = 2(52(y/pl - v^) 2 ) 1/2 p q G M n 



The last expression is related to the Hellinger distance H(p, q) by 
d(p, q) = T*/2H(p, q). The Hellinger distance belongs to the family 
of Cziszar's /-divergences 

D f(p,q) = J f(q/p)dp 

here / is a convex function. As it was shown in pp, the metric given by 
the hessian of /-divergence is a constant multiple of the Fisher metric. 

In the case of a quantum system, the situation becomes more com- 
plicated. In the simplest case, the states of the system are represented 
by density matrices. In analogy with manifolds of classical probability 
densities, a quantum version of the Fisher information metric must be 
decreasing under stochastic maps. Contrary to the classical case, this 
monotonicity condition does not specify the metric uniquely. In fact, 
it was shown by Petz that the monotone metrics can be labelled by 
operator-monotone functions. 

As it was mentioned in 0, there is no general formula for geodesic 
path and distance for a general monotone metrics. Explicit expressions 
are known only in two particular cases, namely the Bures metric and 
the Wigner-Yanase metric. In the present paper, we find an upper 
bound for the geodesic distances for all monotone metrics. This is 
done in a simple way: Following Uhlmann [T5lll9| . we obtain the Bures 
geodesies from certain purifying lifts of curves of density matrices and 
then make use of a duality relation between the smallest (Bures) and 
the largest (RLD) of monotone metrics. It is also shown that this 
upper bound is related to a particular non-commutative version of the 
Hellinger distance. 

2. The manifold and monotone metrics. 

Let M n be the algebra of n by n complex matrices. The set of faithful 
positive linear functionals on M n is identified with the cone of positive 
definite matrices. This set, with the differentiable manifold structure 
inherited from M n , will be denoted by M.. Let D C M. denote the 
sub manifold of density matrices in A4, that is 

V = {peM: Tip = 1} 

The tangent space to M at p G M is T p M = {x G M n : x = x*}. If 
p G V, then the tangent space T p V is the subspace of traceless matrices 
in T P M. 

Let A be a Riemannian metric on M.. Then we will say that A is a 
monotone metric if 



X T(p) (T(h),T(h)) < X p (h,h), P eM, heT p M 
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for all completely positive trace preserving maps T. It is an important 
result of Petz ^H] that a Riemannian metric is monotone if and only if 
it has the form 

X p (h,k) = Tr hJ p {k) 
where J p is given by the operator mean 

(2) J p = R- l [f{L p /R p )]- 1 

Here L p and R p are the left and the right multiplication operator 
and / : (0, oo) — > K is an operator monotone function which is sym- 
metric, that is, fit) = tf{t~ l ). It is immediate from (J2J) that under the 
normalization /(l) = 1, any monotone metric is equal to the Fisher 
metric on commutative submanifolds. Moreover, we have 

2t <m<i±* 



l + t 

for all symmetric normalized operator monotone functions 12 . Ac- 
cordingly, there is a greatest and a smallest element in the set of mono- 
tone metrics. 

The smallest monotone metric is obtained for f(t) = (1 + 1)/2. It is 
called the Bures metric, because it is related to the Bures distance, see 
also Section |U The operator 

J P (h)=g, pg + gp = 2h 

is the symmetric logarithmic derivative, see |§1 IT%1 12]. 

The greatest monotone metric corresponds to the function /(£) = 
2i/(l + t). In this case J p is the right logarithmic derivative (RLD) 

J P (h) = 1 -(p- 1 h + hp- 1 ) 

see [01 EH El • More examples of monotone metrics can be found in 
Section El 



3. Standard representation and monotone metrics. 

The standard representation of the algebra M n is obtained if M n is 
endowed with the Hilbert- Schmidt inner product 

(x, y) = Trx*y 

Let us denote the resulting Hilbert space by H. Then M n is represented 
on H by 

0: M n — > 13(H), a^L a 
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where L a is the left multiplication operator L a w = aw, w G H. Each 
element p in Ai has a vector representative, or purification, w in H, 
such that 

Tr pa = (w, L a w) Va G M n 

Then w G H is a vector representative of p G .M if and only if p = ww*. 

Let pt, t G /, be a smooth curve in . A curve w t in if, such that 
w t is a vector representative of p t for all t G J is called a lift of p t . In 
this case, the tangent vectors are related by 

(3) p t = w t w* + w t Wf 

Let us denote the corresponding projection of the tangent spaces T W H — 
T WW *M by II. 

Let WqWq = p - There are many lifts of p t through w . Among such 
lifts, there is a unique lift with minimal Hilbert space length 

Zff(wt) = J \J (w t ,w t )dt. 

It will be called the horizontal lift. 

The horizontal lift was introduced in [ISJ where the geometric 
phase was extended to mixed states. It was shown that the above 
minimalization problem leads to the condition 

(4) w* t w t = w* t w t 

for all t. The curves wt in H, satisfying this condition, are called 
horizontal curves. The tangent vectors to horizontal curves at w G H 
form a real vector subspace H w = {gw, g = g*}. Let H w be endowed 
with the inner product Re (■, ■), then it is a real Hilbert space, called the 
horizontal subspace. For each h G T WW *M., there is a unique element 
h in H w , satisfying h = U(h). It follows that the inner product in H w 
can be projected onto T WW *M.. As it turns out, this projection defines 
a Riemannian metric on M., moreover 

(5) 4Re(h,k) = 2Trh(L p + R p )- 1 (k), h,keT p M 

is exactly the Bures metric. 

The commutant of <ft(M n ) is the algebra of right multiplication oper- 
ators R a w = wa, a G M n , on H. For each a G Ai, there is an element 
w G H, such that 

(w, R a w) = Tr aa 

This element is given by a = w*w. For each curve a t in Ai, let us 
consider the curves w t in H satisfying w^w t = a t . The tangent vectors 
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of such curves satisfy & t = H(w t ), where IT : T W H — > T W * W M. is given 

II(a;) = x*w + w*x 

We may now proceed exactly as before, choosing for each a t the shortest 
of these curves. It is quite clear that Wt is the shortest curve if and 
only if wl is horizontal, equivalently, w t G H Wt := {w t g, g = g*} for all 
t. Moreover, we have 

(6) x G H w <t=^ x* G H w * 

If we now project the real Hilbert space structure from H w to T W * W A4, 
using the projection II, we will, of course, get the Bures metric again. 
On the other hand, it is easy to see that for each p = ww* and h G T p Ai, 
h := ]^h{w*)~ l is the unique element in H w satisfying h = 11(h). We 
may therefore define 

(7) X p (h, k) := 4Re (h, k) = -Tr p'^hk + kh) 

which is the RLD metric. This shows that there is a duality relation 
between the Bures metric and RLD, see also [HUTU]. 

4. The geodesic distances 

Let A be a Riemannian metric on M. A curve pt, t G [0,1] is a 
geodesic path in M. if its length 

h{pt) = J yj Xptipt, pt)dt 

is the minimum of lengths of all curves connecting p and p\. This 
length is then the geodesic distance of p an d pi- Let us denote by d\ 
the geodesic distance for the metric A in M. and by D\ the geodesic 
distance in T>. 

For Bures metric, the geodesic paths and distances were obtained by 
Uhlmann jTHJ [JJ| as follows. Let po and p\ be two elements in M. and 
let p t be a curve connecting them. If wt is the horizontal lift of pt, then 
by © 

lBuresi.Pt) = 2l H (w t ), 

hence minimizing the Bures length means minimizing the Hilbert space 
length of horizontal lifts of curves connecting po and p\. From the 
definition of horizontality, this minimum is attained at the line segment 
w t = twi + (1 — t)w , such that ||w — Wi\\ is minimal over Wo^o = 
Po, W\w\ = pi. This happens if and only if W\ and w are parallel 
amplitudes, that is, these satisfy Uhlmann's parallelity condition 

(8) w* x w Q > 
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For each u> there is a unique W\ parallel to Wq, given by JH| 

-1/2/ 1/2 1/2x1/2 -1/2 

Wl = Po (Po PiPo ) Po w o 
The geodesic path in Ai, connecting po and p± is then 
p t = {tW! + (1 - t)w ){twi + (1 - t)w )* 
and the geodesic distance is 

d B ures(po, Pi) = 2||iu - u>i|| = 2\Jti p + Tr px - 2Tr (pJ /2 PiPo /2 ) 1/2 

this is called the Bures distance. 

Let now p t be a curve in V, then all lifts of p t are curves on the 
unit sphere S in H. If Wo,^i G 5, the shortest curve connecting them 
lies on the large circle in S through them. The length of such arcs 
for woWq = po and wiw\ = p\ is minimal if wo and w\ are parallel 
amplitudes and, by definition, in this case the arc is also horizontal. 
Hence, the Bures geodesic in D is 

= (w + (1 - t) Wl )(tWQ + (1 - t)wi)* 

Pt \\tw + (i - t) Wl )r 

for parallel amplitudes wo and W\ and the Bures distance 

DBures(Po,Pi) = 2 arccos Tr w wl = 2 arccosTr (pl /2 Pipl /2 ) 1/2 

The duality of the Bures and RLD metrics leads to the following 
upper bound for the RLD geodesic distance. 

Proposition 4.1. Let po, p\ G Ai f then 

— 1/2 n 1/2 

dRLD{Po,Pl) < d B ures(P0,P (Po#Pl) Po ) 

n 1/2/ -1/2 -1/2x1/2 1/2 

Po#Pi = Po (Po PiPo ) Po 
t/ie geometric mean. If po and p\ are in T>, the same holds for 
geodesic distances Drld o,nd D Bures- 

1/2 

Proof. Let wo = p and let w E H be such that wo and w satisfy 
the parallelity condition (jHJ). Then the curve Wt — tw + (1 — t)wo is 
the horizontal lift of the Bures geodesic connecting p and ww*, in 
particular, w t G i?u, t for all t. Then io| is a lift of a curve p t in A^, 
connecting p and w*w and by @, Wj 6 -^™ t *- Consequently, by ((7|), 

dRLD(P0,W*w) < I RLD (Pt) = 2||u>*-U>o|| = 2||w-W || = d Bu res{P0, WW*) 

From the parallelity condition, 10 = gwo for some q = q* > 0. Let us 
choose if such that 

* 1/2 2 1/2 

pi=w w = p > q p ' 
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then q = (p ^ 2 pip 1 ^ 2 ) 1 / 2 and 

* -1/2/ „ n2 -1/2 

WW = p Q (po#pi) p 
The statement for distances in T> is proved exactly the same way. □ 

Remark 4.1. Let Wq, w and q be as in the proof of the previous 
Proposition, then we have 

1 1 wo - w\\ 2 = Trpo + Trpi - 2Trp g 

and 

(9) d Bures (po, Po 1/2 (Pi#Po) 2 Po V2 ) = 2 v /Trpo + Trp 1 -2Trp #pi 
so that po and pi can be exchanged. 

Remark 4.2. Let p* = w* t Wt and g be as in the proof of Proposition 
14.11 Then, in general, pt is not the RLD geodesic. Indeed, it can be 
easily computed that for the RLD metric, the geodesic equation reads 

Pt + T j- = {pt 2 ) - PtPt l pt = a(t)pt 

where a is a smooth function a : / — > R, see also We have 

p t = w*w t = pj /2 (l + t{q - l)) 2 pl /2 

It can be shown by direct computation that the geodesic equation is 
satisfied if and only if 

q(poq - qpo) = (poq - qpo)q 

which, for self-adjoint operators, implies qpo = p$q. It follows that the 
inequality in Proposition 14.11 is strict, unless p and pi commute. In 
that case, the geodesic distances are the same for all monotone metrics. 

In jT2j, a class of generalized relative entropies 

H g (po,pi) = Trpo^Po^PiPo 1 ) 
was introduced, here g is an operator convex function. This is a non- 
commutative version of the /-divergence. It was shown in J7j that the 
generalized entropy H g leads to a constant multiple of the RLD metric 
for infinitesimaly close elements in T>. 

It is easy to see that the right hand side of (JUJ) is equal to ^J2H 90 (pq, pi), 
where 

(10) g {t) = 2 + 2t - At 1 ' 2 . 

Note that on commuting elements, H go is equal to the Hellinger dis- 
tance. 

By maximality of the RLD metric, we get 
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Corollary 4.1. Let po,pi G M. and let X be a monotone metric. Then 



d B ures(po, Pi) < d\(p ,pi) < y 2H go (p , p x ) < 2^/Tr p + Tr p x 
If p ,pi G V, then 

2 arccos Tr {pQ 2 PiPq 2 ) 1 ^ 2 < D\(po,Pi) < 2 arccosTr p #Pi < ^ 

5. The wyd metrics 
The Wigner-Yanase-Dyson (WYD) metrics are defined by 

X a p {h,k) = f a {p + th)U a (p + sk)\ s>t=0 



where 



fa(x) 



1 L ~ a i 1 

-r^X 2 a / 1 

I— a 1 

log(x) a = 1 



As it was shown in jS], these metrics are monotone for a G [—3, 3]. The 
family of WYD metrics is important in quantum information geometry, 
see [TTJ E]. As special cases, for a = ±1, we get the well known 
Bogoljubov-Kubo-Mori metric and for a = ±3 we get the RLD metric. 

The smalest in this family is the Wigner-Yanase (WY) metric, ob- 
tained for a = 0. The WY metric has the form 

X p (h, k) = 4Tr h{ y%, + y%)- 2 (k) 

The corresponding geodesic path and distance was computed in jSj, 
using a non- commutative version of the square root map and a pull- 
back technique. We will show that these can be also easily obtained 
using a similar method as in the Bures case. 

Let p t be a curve in M.. Among its lifts w t w^ = p t , we will again 

choose a horizontal one. In this case, the lift wt is horizontal if it is 

i ii 

contained in the natural positive cone at wo, that is, if Wt = p t uq 
for all t. In this case, the horizontal subspace is = {gu, g = g*}, 
where w = p l l 2 u is the polar decomposition of w. Each tangent vector 
h E T WW *A4 has a unique horizontal lift h° = gu G H®, such that 
h = n(/i°) = gp 1 / 2 + p 1 ^ 2 g- The induced metric 

X p (h, k) = 4Re (h°, k°) = 4Tr h(L l J 2 + R l J 2 ) 2 {k) 

is the WY metric. Note that in this case x G i?° if and only if x* G if ° * , 
so that the WY metric is self-dual, in the sense mentioned in Section 
01 Let us also remark that it is possible to obtain all the monotone 
metrics in a similar manner, see |H QUI . 

Let now po and p\ be in M. and let p t be a curve connecting them. 
Again, the WY lenght of p t is twice the Hilbert space length of its 



1/2 

horizontal lift u>t = Pt u o- Therefore, pt is the geodesic path if wt 
tpl^ito + (1 — t)p^ 2 Uo, that is 

Pt = (tp\ /2 + (i-t)pl /2 f 

and the geodesic distance is 



<W(po,Pi) = 2||pJ /2 - pl /2 || = 2^/Trpo + Tr Pl - 2Ti p] 12 p\ /2 

Similarly, if po,Pi G V, then p t is a geodesic path if and only if w t lies 
on the large circle connecting p l J 2 u and pj^Mo- Hence 

(tpl /2 + (l-t) P l /2 ) 2 
Pt \\tp\ /2 + {l-t)pl /2 \? 

and 

D W y(Po,Pi) = 2arccosTrpJ /2 p} /2 

Let us denote by A^ p = L^R^ 1 the relative modular operator. In 
[To] , a class of quasi-entropies was introduced by 

S g (p,a) = Trp^g(A^)(p^) 

where g is an operator convex function. This is another quantum ver- 
sion of the /-divergences. It is easy to see that 



dwy(Po,Pi) = \j2S go (po,pi), 



where go is given by (fTT)|) . It was proved in ^Hj that each monotone 
metric can be obtained as the hessian of S g for a suitable operator 
convex function g. The choice g = go leads to the WY metric. 

From the previous Section and the fact that the WY metric is the 
least element in the family of WYD metrics, we get 

Corollary 5.1. Let X be a WYD metric and po, Pi G Then 



^2S go {po,pi) = d WY {p , Pi) < d x (p , pi) < ^2H go (po,Pi) 
where g (t) = 2 + 2t - At 1 / 2 . If p , pi G V, then 

2 arccos Tr p^ 2 P\ 2 < -Da(PO)Pi) < 2 arccos Tr p #Pi 
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